Wandering breathers and self-trapping in weakly coupled nonlinear chains: classical 
counterpart of macroscopic tunneling quantum dynamics 
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We present analytical and numerical studies of phase-coherent dynamics of intrinsically localized 
excitations (breathers) in a system of two weakly coupled nonlinear oscillator chains. We show that 
there are two qualitatively different dynamical regimes of the coupled breathers, either immovable or 
slowly-moving: the periodic transverse translation (wandering) of low-amplitude breather between 
the chains, and the one-chain-localization of high-amplitude breather. These two modes of coupled 
nonlinear excitations, which involve large number of anharmonic oscillators, can be mapped onto 
two solutions of a single pendulum equation, detached by a separatrix mode. We also show that 
these two regimes of coupled phase-coherent breathers are similar, and are described by a similar 
pair of equations, to the two regimes in nonlinear tunneling dynamics of two weakly linked Bose- 
Einstein condensates. On the basis of this profound analogy, we predict a new tunneling mode 
■ of two weakly coupled Bose-Einstein condensates in which their relative phase oscillates around -| 

modulo 7r. We also study two-chain breathers, which can be considered as bound states of discrete 
breathers with different symmetry and center locations in the coupled chains, and bifurcation of 
the anti-phase two-chain breather into the one-chain one. We show that static displacements of 
the coupled chains with nonlinear localized excitation, induced by cubic term in the intra-chain 
anharmonic potential, scale approximately as the total vibration energy of the excitation, either 
one- or two-chain one, and do not depend on the inter-chain coupling. Delocalizing transition of 
ID breather in 2D system of a large number of parallel coupled nonlinear chains is described, in 
' which the breather, initially excited in a given chain, abruptly spreads its vibration energy in the 

whole 2D system upon decreasing breather frequency or amplitude below the threshold one. The 
threshold breather frequency is above the cut off phonon frequency in 2D system, and the threshold 
breather amplitude scales as square root of the inter-chain coupling constant. Delocalizing transition 
of discrete vibrational breather in 2D and 3D systems of coupled nonlinear chains has an analogy 
with delocalizing transition for Bose-Einstein condensates in 2D and 3D optical lattices. 

> 
in 

Nonlinear excitations (solitons, kink-solitons, intrinsically localized modes and breathers) can be created most easily 
in low-dimensional (ID and quasi-lD) systems [l|, 0, 0, 0, 0, 0, 0, 0, 0, OH- Recent experiments have demonstrated 
the existence of intrinsically localized modes and breathers in various systems such as 2D and 3D arrays of nonlinear 
r*"** . optical waveguides [TJ Elj], low-dimensional crystals [l3|], antiferromagnetic materials [HI], Josephson junction arrays 
[l5| . photonic structures and micromechanical systems [l7j . protein a- helices [l8| . and a-uranium D3- Slowly- 
moving breathers and supersonic kink-solitons were also described in ID nonlinear chains [is mm mm, dna 

One-dimensional arrays of magnetic or optical microtraps for Bose-Einstein condensates (BECs) of ultracold quan- 
tum gases with tunneling coupling provide a new field for the studies of coherent nonlinear dynamics in low-dimensional 
systems [13, In the mean- field th eory , the tunneling coupling between two BECs is similar to the linear cou- 
pling between two optical waveguides [29| or between two oscillator chains (nonlinear phononic waveguides). Here 
we show that phase-coherent dynamics of macroscopic ensembles of classical particles (weakly-localized breathers) in 
two weakly linked nonlinear oscillator chains has a profound analogy, and is described by a similar in every respect 
pair of equations, to the tunneling quantum dynamics of two weakly linked BECs in a macroscopic double- well po- 
tential (single bosonic Josephson junction) (3(|. The exchange of energy and excitations between the coupled classical 
oscillator chains takes on the role which the exchange of atoms via quantum tunneling plays in the case of coupled 
BECs. Therefore such phase-coherent energy and excitation exchange can be considered as a classical counterpart of 
macroscopic tunneling quantum dynamics. 



I. INTRODUCTION 
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We show that there are two qualitatively different dynamical regimes of the coupled breathers, the oscillatory ex- 
change of the low-amplitude breather between the chains (wandering breather) , and one-chain-localization (nonlinear 
self-trapping) of the high-amplitude breather. These two regimes, which are detached by a separatrix mode with zero 
rate of energy and excitation exchange, are analogous to the two regimes in nonlinear dynamics of macroscopic con- 
densates in a single bosonic Josephson junction [30(. Essentially, the phase-coherent dynamics of the coupled classical 
breathers is described by a pair of equations, which coincides with the pair of coupled mean-field equations describing 
coherent atomic tunneling in a single bosonic tunnel junction [3lL |32|. The predicted evolution of the relative phase 
of two weakly coupled coherent breathers in both dynamical regimes is also analogous to the evolution of relative 
quantum mechanical phase between two weakly coupled macroscopic condensates, which was directly measured in a 
single bosonic Josephson junction by means of interference [30]. Moreover, the predicted separatrix in the excitation 
exchange between macroscopic ensembles of phase-coherent particles (weakly localized breathers) in coupled oscilla- 
tor chains with "repulsive" nonlinearity, which is determined by the ratio of the intra-chain nonlinearity (intra-chain 
interaction) and the inter-chain coupling, can be considered as a classical nonlinear dynamical model of the reversible 
interaction-induced superfluid-Mott-insulator quantum phase transition, observed in BEC with repulsive inter-atomic 
interaction in a lattice with tunneling inter-site coupling [33]. All these results bring to light a striking similarity, 
both in display and evolution equations, between the classical phase-coherent excitation exchange and macroscopic 
tunneling quantum dynamics which can motivate new predictions and experiments in both fields. On the basis of 
this profound analogy, we predict a new tunneling regime of two weakly coupled BECs in which their relative phase 
oscillates around ^ modulo tt, which can be observed by means of interference. This regime is different from the 
regime of Josephson plasma oscillations, already realized in experiments (30j , when the relative phase of two weakly 
linked BECs oscillates (or fluctuates [3J]) around zero. 

We also show that static displacements of the coupled chains with nonlinear localized excitation, induced by cubic 
term in the intra-chain anharmonic potential, scale approximately as the total vibration energy of the excitation, either 
one- or two-chain one, and do not depend on the inter-chain coupling. This observation is also valid for a narrow stripe 
of several parallel coupled nonlinear chains. We also study two-chain breathers which can be considered as bound 
states of discrete breathers with different symmetry and center locations in the coupled chains, and bifurcation of the 
anti-phase two-chain breather into the one-chain one. Bound states of two breathers with different commensurate 
frequencies are also found in the two-chain system. Merging of two breathers with different frequencies in one breather 
in two coupled chains is described. Periodic transverse translation (wandering) of low-amplitude breather in a system 
of several, up to five, coupled nonlinear chains is observed, and the dependence of the wandering period on the number 
of chains is analytically estimated and compared with numerical results. Delocalizing transition of ID breather in 2D 
system of a large number of parallel coupled nonlinear chains is described, in which the breather, initially excited in 
a given ID chain, abruptly spreads its vibration energy in the whole 2D system upon decreasing breather frequency 
or amplitude below the threshold one. The threshold breather frequency is above the cut off phonon frequency in 
2D system, the threshold breather amplitude scales as square root of the inter-chain coupling constant, and breather 
vibration energy is localized mainly in one chain at the derealization threshold. Similar delocalizing transition for 
ID breather should occur also in 3D array of parallel coupled nonlinear chains. Delocalizing transition of discrete 
breather in 2D and 3D systems of coupled nonlinear chains has an analogy with delocalizing transition for polarons 
in 2D and 3D lattices [35| and for Bose-Einstein condensates in 2D and 3D optical lattices [36[ . 




The paper is organized as follows. In the next Section II we describe our model and the analytical predictions, 
derived on the basis of this model. In Section III we describe our numerical simulations and comparison with the 
analytical predictions, which include simulation of dynamics of discrete breathers in two weakly coupled chains (Section 
III A), simulation of wandering breather in two weakly coupled chains (Section III B), and simulation of breathers in 
a system of M > 2 parallel weakly coupled anharmonic chains (Section III C). In Section IV we give a Summary of 
the main new results in this paper. 



We consider two linearly coupled nonlinear oscillator chains (with unit lattice period). We model the chains with 
the Fermi-Pasta- Ulam (FPU) Hamiltonian, which is one of the most simple and universal models of nonlinear lattices 
and which can be applied to a diverse range of physical problems [37j : 




II. MODEL AND ANALYTICAL PREDICTIONS 
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Here is displacement of the n-th particle from its equilibrium position in the i-th chain, p$ = v!$ is particle 
momentum, jW, cS l \ fiW and C are, respectively, intra-chain linear, intra-chain nonlinear and inter-chain linear force 
constants. We assume that the coupling is weak, C^l"> , and do not include the nonlinear inter-chain interaction. 
The /3-FPU Hamiltonian |T]) (with aft = 0) describes, in particular, purely transverse particle motion Q. Torsion 
dynamics of DNA double helix can also be approximated by the /3-FPU Hamiltonian (1) [24[. On the other hand, 
weakly coupled nonlinear molecular chains in polymers are characterized by the asymmetric intra-chain anharmonic 
potential, with nonzero [25L f26j] . 

Hamiltonian (1) generates the following equations of motion, z = 1,2: 



u 



Below we consider two chains with identical anharmonic force constants, when = = a and (3^ = (3^ = 
(3. We are interested in high-frequency and therefore short-wavelength dynamics of the coupled chains, when the 
displacements of the nearest-neighbor particles are mainly anti-phase. For this case we can introduce continuous 
envelope-functions for the particle displacements in the chains, fn^ = u„\— 1)", f^ 1 = f(x)i- These envelope- 
functions f(x)i are supposed to be slowly varying on the interatomic scale in both chains, dfijdx « 1, which allows 
us to write corresponding partial differential equations for these functions, see, e.g., Refs. 0, @, 0, Hy. To obtain 
nonlinear envelope-function equations for the coupled chains with asymmetric inter-particle potential, one needs to 
consider separately the dynamics of the relative, = u^ +1 — u„ , and total, Wn = u n+i + u n\ displacements of 

the nearest-neighbor particles, cf. Ref. 2]. Equations for v$ and w$ can be obtained from Eqs. ^) and look as 
following: 

= l ( nv«li+vV 1 -2vU 1 ) + a[vW +v M-2vW] 

+ - 2v^] + C(v^ - „«), (3) 

+P[vW-vV 3 1 ] + C(w£-»-wV). (4) 
Introducing continuous relative and total displacements, w„ = v(x)i and u>„ = w(x)i, and making in Eqs. ([3]) 

. . u\ u\ U~\ (£\ 

and (|4J) expansions of the differences u n ± 1 — u„ and w„±i ~ w n up to the second order, we get the following partial 
differential equations for v(x,t)i and w(x, t),: 

fi < = -l^^-^ ) Vi-^(av i +/3v?)-4f3vf + C(« 3 -i-^), (5) 

^ = l W^ + 2 a^+2^ + C(w 3 - i -w i ). (6) 

The main small parameters, which will allow us to find asymptotic solutions of these equations, are the amplitudes 
of the displacements Mn' of the coupled oscillations and hence their envelopes f(x,t){. As follows from Eqs. ((6|), in 
the chains with asymmetric inter-particle potential (with finite a) the total displacements Wi of the nearest neighbors 
are weakly coupled (via anharmonic interaction) with relative displacements Vi and are in general small in comparison 
with them, Wi « Vi (see Eqs. (9)-(ll) below). In such chains, the asymmetric potential induces both the static 
(zero-harmonic) and second-harmonic terms in the total displacements of the nearest- neighbor particles, but the latter 
terms are smaller than the former ones in the small-amplitude limit, cf. Ref. [2]. 

Now we analyze the static displacements woi(x) of two identical weakly coupled chains with = l^ = 1. In 
order to minimize the static elastic energy, the static displacements of the chains should be equal at both infinites: 
wqi(x) — Wq2(x) for x — > — oo and x — > +oo. To find WQi(x), we take the sum and the difference of Eqs. © for the 
two chains to obtain the following equations: 

d 2 { Wl +w 2 ) „ d(vl+vp 0R d(vf + vl) 
OX ox ox 

= - 2c iwi - W2) + 2* d Mpt + 20 Mp31. (8) 

ox ox ox 

Omitting the time-derivative terms, we get from Eq. ((7]) the following expression for the static center-of-mass 
displacement of the two chains: 

w( cm ){x) = %m^m = a r {vl + vl)d ^ (9) 



m + W2 = - ^ ; W2J + 2a ^ Ul Z v *' + 2p " yUl Z V2 > , (7) 
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where the angle brackets denote the time averaging. The nonlinear term, proportional to (3{v\ + uf), cf. Eq. (J7J), is 
averaged out in Eq. (9) because of the third power of yW. At both infinities, the displacement w\^ M \x) coincides 
with equal static displacements of the coupled chains, w\^ M ^ (±00) = '{±00) — Wq 2 ^(±oo). It is assumed in 

Eq. (9) that Wq (—00) = 0, v^^—oo) — and W q x — (—00) = 0, i = 1,2. Therefore the static displacements 
of the both coupled chains Wq(x)^ have in general the form of the kink-like pattern with equal center-of-mass 
displacement difference Awq CM ' = Wq CRI \+oo) — Wq CM \—oo) = Wq CM \+oo). The kink-like static displacement 
pattern of intrinsically localized mode (breather) in a single a-/3-FPU chain was studied in Refs. (2I l38l. l39l. l40l| . 

The term in angle brackets in r.h.s. of Eq. (9) coincides with the total vibration energy of two coupled chains (in 
the low-amplitude limit when the energy is governed mainly by the quadratic terms in Hamiltonian (1)). Therefore 
the center-of-mass displacement of the coupled chains with asymmetric anharmonic potential scales approximately 
as the total vibration energy of the nonlinear localized excitation, the one- or two-chain breather in two coupled 
chains: w^ M ^ — —a f^°°{v 2 + v 2 )dx . Equations (7) an (9) also show that static strain of the chains, given by the 
displacement derivative, is determined by the local density of the total energy of the nonlinear localized excitation: 
dwf M \x)/dx = —a(vf + v|)(a;). These conclusions are confirmed by our numerical simulations, see Figs. 2(c) and 
2(d) below. 

These conclusions can be easily generalized to the system of M > 2 parallel coupled nonlinear chains, see also 
Section III C below. In this case only the last term in Eqs. (5) and (6), describing the inter-chain coupling, will be 
modified (cf. Eq. (62) below). In a narrow stripe of chains, when the number of particles N in each chain greatly 
exceeds the number of chains, N ^> M , the static center-of-mass displacement in the multi-chain system will have 
the form: 



1 m icv r x 



M „ „ x M 

(CM) 

M ^ -<"»v"v M 

where v m (x) and WQ rn (x) are, respectively, the relative and static total nearest-neighbor displacements in the m-th 

chain. It is assumed in Eq. (10) that Wq (—00) = 0, v m {— 00) = and w ° dx — (—00) = 0. As follows from Eq. (10), 

the lattice displacement w ( ^ M \x) and the corresponding lattice strain dw^ 1,1 ^ jdx in the multi-chain system depend 
on the number of the excited chains and their excitation energies. Importantly, the static lattice displacements and 
strain in a system of coupled chains, Eqs. (9) and (10), do not depend on the strength of the inter-chain coupling, 
which is connected with long-range action of static elastic fields. 

The magnitudes of the static displacements in the central region of the kinks in the two chains are different in 
general and can be found with the help of Eq. (8) and the above requirements for the values of wo(x)^ at the 
infinities. Omitting the time-derivative terms, we get from Eq. ([8]) the following solution for the difference of static 
displacements of the two coupled chains: 

1 /Y f x 

Aw (x) = w i(x) - wq 2 (x) = -W-( / sinh[(s - x)V2C]A(s)ds + 

4 V C J _ co 

+ / sinh[(a; - s)V2C]A(s)ds - / sinh[(s + x)V2C]A(s)ds + 

(11) 



/OO 
smh[(s + x)V2C]A{s)ds). 
-X 



where A(s) — ad(v1(s) — V2{s))/ds. Taking into account that vf(x) is a symmetric function with respect to the 
breather center x — 0, see Eqs. (20) and (21) below, Eq. (11) gives us that Awo{x) — -Aw (-i) and therefore 
Awo(x) is zero both at x = and at both infinities. [As follows from Eq. (11), the difference of static displacements 
of the coupled chains Awq(x) depends in general on the coupling constant C, in contrast to the static center-of- 
mass displacement Wq CM \x), Eq. (9).] As follows in turn from Eqs. (6), the second-harmonic terms in the total 
displacements w{x) 1 in the chains have the smallness of the third order (in v(x)^) because of the smallness of spatial 
derivative dv^ 2 /dx ~ Xiv^ 2 ~ v^ 3 , see Eqs. (19), (20) and (23) below. This allows us to neglect in the following 
analysis of breather dynamics in coupled chains both the difference of static displacements of the chains Awq (x) and 
the second-harmonic terms in the total displacements w(xY in comparison with the static center-of-mass displacement 

wf M \x).^ 

Taking into account that v{x)i — 2f(x)i, we get from Eqs. (5), (7) and (9) the following coupled nonlinear partial 
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differential equations for the real envelope-functions f(x)i, i = 1,2: 



o2 t 

I + u? mi U + —^ + 16/3/f - 16a 2 (A 2 + /I)/, - C/g-i = 0, (12) 



where u) m i — y/Al^ + C is characteristic frequency slightly above the maximal phonon frequency in the i-the isolated 
chain. These self-consistent equations describe the dynamical relative lattice displacements Vi(x), while the static 
lattice displacements woi(x) can be obtained with the use of Eqs. (9) and (11). It is worth mentioning that collective 
static deformation of the coupled chains, induced by the asymmetric intra-chain anharmonic potential, makes the 
coupled envelope- function equations (12) non-local in general, via the term proportional to —a 2 . 

In order to deal with the amplitude and phase of the coupled nonlinear excitations, it is useful to introduce complex 
wave fields ^(x, t)i for each chain, see, e.g., Ref. [2l| : 

f(x,t) i = ±[*(x,t) i + *(x,t)*]. (13) 

Then in resonance approximation we get from Eqs. (12) and (13) the following coupled equations for ^(x,t)i, 
i = 1,2: 

+ + <*&i*<) + W I *< I' *i ~ 8 « 2 d *i I' + I *2 | 2 )*< - ^C# 3 -i - 0, (14) 

(and complex-conjugated equations for ^*). This approximation, in which we neglect the higher harmonics, is valid 
for the dispersive system under consideration due to the weakness of nonresonant interaction between the mode with 
fundamental frequency and its third harmonic, cf. Ref. [2ll |. 

Using Eqs. (14) in the assumption that the shift of the wandering breather frequency, caused by nonlinearity, weak 
inter-chain coupling and slow breather motion along the chains, is relatively small, one can readily show the existence 
in general of the following three integrals of motion: 



N b = J [| *i | 2 + | * 2 | 2 ] da, (15) 

2 

E b = I I *< I 4 -8a 2 (| * x | 2 + | * 2 | 2 ) | * 4 | 2 - 

2—1 

1 1 

~2 ) ~2 C( * 1 * 2 ^ VM/ * )]dx ' (16) 



Pbx = 4/^ 



1 \ T *d*i 

— « — 

dx dx 



dx, (17) 



which play the role of breather total number of excitations, total energy and total momentum along the chain axis, 
respectively, cf. [2l[. The existence of these integrals of motion demonstrates that the exchange of energy and 
momentum between two nonlinear systems (1) is a coherent phenomenon, which depends in general on the initial 
excitation conditions. 

It is also useful to define partial numbers of excitations 

Ni = J |*i \ 2 dx, »=1,2, 

when Ni + Ji — 0, where 

Ji = -Ja-i = ^ / [***3-< - *?*3-<] dx (18) 

is a total inter-chain flux of excitations, which conserves the total number of them: Nb = where Nb = N± + N2, u 
is frequency of wandering breather, see Eqs. (19), (20) and (22) below. 

To describe with the help of Eqs. (14) a slowly-moving breather, wandering between two weakly coupled nonlinear 
chains with positive (repulsive) anharmonic force constant f3, we assume the following form for the complex fields ^1 
and "J^: 

exp\i(kx — tut)] _ . i . „. 

*i = * max ttt , ^cosGexp (--$), 19 

cosh[Ai(a; — Vt)\ 2 

exp\i(kx — u)t)] „ , i . , 

*2 = *ma X * , ( ^sinGexp (-$ , 20 

cosh[A2(2; — Vt)\ 2 
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where uj, V <C 1 and k <C 1 are, respectively, frequency, slow velocity and small wavenumber related with the moving 
breather, A1.2 <C 1 are real parameters describing inverse localization lengths of the breathers; $ = $(t— kx/uj) stands 
for the relative phase of the coupled chains, while the parameter O = Q(t— kx/ui) describes the "relative population" 
of the two chains z = (n\ — ri2)/(ni +TI2) = cos 20, where nj=|\l/i| 2 is local density of excitations in the i-th chain, and 
(| ^1 | 2 + I *2 | 2 ) = ^max m Eq. (14). The parameters 6 and $ determine also the inter-chain flux of excitations, 
cf. Eq. (18): 

T j = C^max f sin 29 sin $ 

1 2 2uj J coBh[/si(a:-Vt)]oosh[K2(a;-yt)] 1 ' 

Using Eqs. (14), (19) and (20), after some algebra we obtain dispersion equations for the introduced parameters, 

^ = l("li+^2) + (W~8a 2 )^ ax -k 2 -c£^- y (22) 

\\ = (G/We-Sa 2 )* 2 ^, \ 2 2 = (Wsm 2 G-8a 2 )* 2 max , (23) 
duo , . 

= m> (24) 

and evolution equations for <E> and O: 

4> = ^(c, 2 ll -^ 2 ) + ^™^cos(2e) + -cos$cot(2e), (25) 

LUJ UJ UJ 

C 

= — sin$. (26) 
2uj 

In the derivation of Eqs. (j2"5|) and (j2"6")l . it was assumed explicitly that the ratio is equal to one. 

The latter is valid for small-amplitude breathers with long localization lengths, A12 -C 1. In this case the above 
assumption, which is exact for the central region of the breathers, x-Vt& 0, will be (approximately) valid for a large 
number of particles, which form weakly localized wandering breather in weakly coupled nonlinear chains. It is also 
assumed in the considered approximation that the shifts of breather frequency uj, Eq. (22), caused by weak inter-chain 
coupling C, nonlinearity (3/3 — 8a 2 )^ nax and slow breather motion along the chains fc 2 , as well as the characteristic 
frequency difference \uj m i -uJm2\/uJ, are all relatively small. It is worth to mention that the parameter — 8a 2 5' 2 la:I . 
determines the shift of breather frequency uj and inverse localization lengths Ai,2, see Eqs. (22) and (23), but does not 
enter evolution equations (25) and (26) for and $. Equations similar to Eqs. ([2!>]) and (j2"B")) were probably derived 
for the first time in Ref. [41( for the description of energy exchange between two weakly coupled classical anharmonic 
oscillators (with uj m \ — uj m 2 and a = 0). Later similar equations were used for the description of the power exchange 
between two weakly coupled nonlinear optical waveguides [42| and of the dynamics of coupled nonlinear oscillators 

63- 

Equations (25) and (26) can be written in an equivalent form for the relative phase $ and population imbalance z, 
when z — cos 20 and \/l — z 2 = sin 20: 

$ = uj ml - uj rn2 + 3(3 *™* z + - 2 cos $, (27) 
UJ uj VI - z 2 

i = \J\ - z 2 sin$. (28) 

UJ 

Here the variables $ and z are canonically conjugate, 

• dH e ff . dH e ff , 

* = ^f' Z = ~^f> (29) 

with the following effective Hamiltonian (which has the dimension of frequency) : 

= 3^L, Z 2 _ C y / 1 _ z 2 cosq> + z{LUml _ Wm2 ). (30) 

ZUJ UJ 



The very same equations (|27|) and (|28|) for <I> and z, which are equivalent to Eqs. (|25[) and ([26]) for $ and 0, were 
obtained in [3l], [13] m the mean-field theory of tunnelling dynamics of two weakly coupled Bose-Einstein condensates, 
which were later used in the analysis of the experimental realization of a single bosonic Josephson junction [30| . 
Therefore generic evolution equations (27) and (28) (being written in corresponding units) for the dynamics of two 
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coupled phase-coherent breathers or tunneling dynamics of two coupled BECs do not explicitly depend on the " source" 
equations: Eqs. (14) for coupled oscillator chains or mean-field Gross-Pitaevskii equations for coupled BECs (3ll.[3^|. 
In our case, Eqs. (|27| and ([28]) describe the exchange of lattice excitations between the chains rather than atomic 
tunneling. One can therefore consider such excitation exchange as a classical counterpart of tunneling quantum 
dynamics. Some other equivalent forms of the effective Hamiltonian (30) were also discussed in Refs. [44 145|. 

It is noteworthy that equations, similar to Eqs. ([25]) and (126|) . describe the dynamics of two weakly coupled nonlinear 
oscillators (with different harmonic eigenfrequencies ui m i and io m 2) [43l l46l|. Therefore wandering breathers can be 
considered as weakly coupled phase-coherent nonlinear macroscopic oscillators. 

Equations (25) and (26) can be solved analytically for the given initial conditions with the use of variational method. 
For two identical chains, with cj ml =cj m2 =u; m , we assume that the solution of Eq. (26) has the following form: 

cos$ = A(t)/sm(26), (31) 

where A=0 for sin(29)=0. Using this ansatz, we get from Eqs. (25) and (26) that 

A = -G ^™* sin(26) cos(26)6. (32) 
For the aforementioned initial condition, when A=0 for sin(20)=O, from Eq. (32) we find: 

A = -^^sin 2 (26), (33) 

2 Q/?iTf2 



cos$ = - 3 ^ ax sin(29) = - 3 ^" lax V±-z 2 , (34) 
which corresponds to <f> = 7r/2 modulo 7r for sin(29) = 0. This exact solution of Eqs. (J35J) and (J2SJ), $27} and (gSJ), 

33 *]/ 2 

conserves the effective Hamiltonian: H e ff— 2 ™ ax ■ The important feature of this solution is that the relative phase 
$ is self -locked to the value ^ modulo 7r by the total population imbalance | z |=1 of the two coupled chains. The 
phase portrait of Eq. ([34]) in the $ — z plane is given by 

(kcos$) 2 +z 2 = 1, (35) 

where k = 2C/3/3* 2 naa; , see Fig. 1. 

With the use of Eq. {34]), Eqs. (J22j> , ([25]) and (J26]), J37J) and (H) will take the following form: 



^ = co 2 m + (^(3-8a 2 )^ 2 nax -k 2 , (36) 

^ = 1^/3*^08(26), 6=^sin<I>, (37) 
3 „„, C 



* = ir^naxz, z = VT^sin*. (38) 

lui to 

Finally, for two identical weakly coupled chains with l^ = l^> = 1 and io m ~ 2 we get from Eqs. (37) the following 
two equivalent physical-pendulum equations: 



for 5 = 46 where f2 = 3/3* 2 naa ,/4, and 



E + fi§sinE = (39) 



C 2 

8 + — sin 5 = (40) 
4 



- for S = 2$ — 7r. In the following we will solve Eq. (139[) with the initial condition 6(0) = 0, which corresponds to 
zero complex field ^2 in the second chain at t=0 (or z(0)=l), and which is realized in our simulations. Therefore we 
will assume that S(0) = 0, $(0) = ir/2 modulo it and 2(0) = C. The corresponding initial conditions for 5 in Eq. 
(40) are 5(0) = modulo 2tt and 6(0) = \{i^ 2 max . 

Solution of the nonlinear physical-pendulum equation is well known, see, e.g., Ref. [471 ]. and can be written in terms 
of elliptic functions (with the elliptic modulus k=2C /3(3fy 2 nax ) . Namely, one can obtain from Eq. (39) the following 
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time evolution and fundamental oscillation frequency uj^q for 5 = 49: 

49 = 2 arcsin [Ksn(fV, «)], « < 1, (41) 
49 = 2 arcsin[s7i(Kfi *, !/«)]> « > 1, (42) 

" 4e = 2^00' 16 - 1 ' (43) 

uj4& = n Km' ^ (44) 

where sn is the Jacobian elliptic sine, K(k) = F(it/2,k) is complete elliptic integral of the first kind. The solution 
for 6 = 2$ — 7r can also be found with the use of Eq. (40) . 

The important property of the solution of the physical-pendulum equation (39) (or (40)) is the existence of two 
qualitatively different dynamical regimes of Eqs. (37) and (38), which are detached by a separatrix corresponding to 
the condition n = 1 or 3(3^^^ = 2C. 

For k3>1 or 0^ nax '^.2C/3, the parameter 9 linearly grows with the "running" time t=t— while <& « $(0) = 5 
(see, e.g., Ref. gjj): 

8 » C- + p J" aa sin(Ct), ™ ax sin(C-), 2 » cos(C-), C = C--^ = C- P ' naa . (45) 

4 64C 2 2 2C 2 ; ' v 2 ; ' C 16C v 7 

In this regime 9 spans the full range from to 2n, which, according to Eqs. (19) and (20), corresponds to the 
complete energy exchange between the nonlinear chains and therefore to the breather, wandering between the two 
chains. According to Eq. (21), the flux of the inter-chain energy exchange in this regime J\ oc C^ max sin 29 sin A = 
C^ m ax sin(Ct/2) has the rate, u>beat = C/2 = (C — 9/3 2 5 ,4 aaa ./(16C'))/2 for k^>1, which continuously decreases with 
the increase of the ratio P^ 2 max IC below the separatrix. Since (cos9 2 ) = ^sin9 2 ) = 1/2 in this mode, the time- 
averaged inverse localization lengths A1.2, Eq. (23), are equal in the two chains: Ai = A2 = ^maxy/SP — 8a 2 (for 
3/3 > 8a 2 , see Section III A below). 

This dynamical regime is analogous to the complete energy exchange in the beating in a system of two coupled 
harmonic oscillators. The solution 9 = Ct/4, $ = ■§ and z = cos(Ct/2) can be obtained directly from the linearized 
Eqs. (14). 

Similar tunneling mode can also be realized for BEC in a symmetric double-well potential when the condensate is 
initially loaded into one of the wells, z(0)=±l, cf. Ref. [29). In such mode the relative phase of the coupled BECs will 
oscillate around \ modulo tt, see Eq. ()45|) and Fig.l, which can be observed by means of interference. This tunneling 
mode is similar to but is different from the mode of Josephson plasma oscillations, already realized in experiments 
[30l ] . in which the relative population imbalance is always less than one and the relative phase of the coupled BECs 
oscillates around zero (modulo 2tt). The observable difference between the time evolution of average relative phases 
in the observed in 30] tunneling mode and the mode predicted in Fig. 1 is caused by the difference in the initial 
states (population imbalance and relative phase) of two weakly linked Bose-Einstein condensates: it is | z(0) \< 1 and 
$(0) = modulo 2ir in the former case while it is | z(Q) |= 1 and $(0) = ir/2 modulo ir in the latter case. 

It is worth to emphasize that according to Eqs. (19), (20) and (45), the temporal Fourier spectrum of the wandering 
breather is determined both by the frequency uj and time dependence of 9. The latter produces in spectral density 
of the wandering breather a series of doublets of frequencies, shifted upward and downward with respect to the 
fundamental breather frequency ui: lo ± nw(, eat = to ± nC/2, n = 1,2, see Eq. (45). The lower frequencies in such 
doublets can become less than the threshold frequency u> rn in the phonon band of the coupled chains, see the dashed 
line in Fig. 8 below, which results in weak damping of the coupled nonlinear oscillations caused by phonon emission 
in the chains. 

In the opposite limit fc C 1 or P^max ^> 2(7/3, one has 

9 w — 2- — sin(7i3$L T t), $ « ~ + ^-P^ 2 n J H 2 sin(^/3^ 2 n J), 

3P^ 2 4 2 4 18(/3\E f2 ) 2 2 

2C 2 3 C 2 

■ S m 2 (-p* 2 max t), ^ 2 max = ^ 2 max -——. (46) 



9(^v|/2 \2 v 4 max " ^ max ^ max 9/M 2 



max 



In this mode the inter-chain energy exchange between the coupled chains is incomplete since always 9 <C 1, cf. 
Eqs. (19) and (20). With the increase of the ratio P^ max /C beyond the separatrix, the flux of such incomplete 
energy exchange, J\ oc (C 2 /^ ma2: ) sin29 sin $ = (C 2 /2^ max ) sin(3p^l 2 nax t/2), cf. Eq. (21), gradually decreases but 
its rate 3P^ 2 nax /2 gradually increases, see Eq. (46). This mode is similar to the macroscopic quantum self-trapping 
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of Bose-Einstein condensate in a single bosonic Josephson junction [2 



3( 



3ll |32{ |. as well as to the asymmetric 



nonlinear mode (known, e.g., for two coupled nonlinear waveguides |42T 4917 5C | ) . in which one system, here is system 
1, carries almost all vibrational energy while the other one is almost at rest. 

The separatrix k = 1 or ft^ 2 max — 2C/3 is characterized by zero frequency of the physical pendulum (39) (or (40)), 
which corresponds to the infinite period of inter-chain energy exchange. For the considered initial conditions, the 
separatrix is described by the following solution of Eqs. (37) and (38): 

C ir C 7r 

9 = arctan[exp(— t)} - -, $ = 2 arctan[exp(— t)], z(0) = 1, 2(00) = 0, $(0) = -, $(00) = n. (47) 

The flux of the inter-chain excitation exchange is zero at the separatrix for t — ► 00: Ji ~ since sin$ ~ 0, see Eq. 
(21). On the other hand, the breather frequency uj is finite at the separatrix and is equal to 

u> sep = u m + (9/3/8 - 2a 2 )* 2 max - fc 2 /4 « 2 + C - 2a 2 * 2 nax - k 2 /A, (48) 

since 3(3^^^ = 2C at the separatrix and ui m w 2 + C/4, see Eqs. (12), (19) and (36). The value uj sep ~ 2.1 for 
C = 0.1, a — and k <C 1 coincides with good accuracy with the one, oj sep ~ 2.098, which follows from our numerical 
simulations of the separatrix-like dynamics of the wandering breather in two coupled /3-FPU chains, see Figs. 7 and 
11 below. In corresponding dimensionless units, the above separatrix, Eq. (47), gives the position, in terms of spatial 
soliton amplitude and inter-fiber coupling, of the numerically revealed separatrix-like regime between the regimes 
of the total and partial exchange of photon energy between two coupled nonlinear optical fibers (fiber directional 
coupler), described by two coupled nonlinear Schrodinger equations [5l|. The existence of a separatrix was revealed 
in different models of nonlinear dimer (two coupled nonlinear systems with two degrees of freedom), both classical 
and quantum, in Refs. 0, [H, HH, H3, [53| . 

The important result of our studies in this field is that the dynamics of two coupled nonlinear systems with many 
degrees of freedom (two weakly coupled low-amplitude and macroscopically wide breathers) can be (approximately) 
mapped onto the dynamics of exactly solvable nonlinear dimer. The (macroscopically) large number of degrees of 
freedom in the considered coupled systems also allows one to consider their relative phase <I> as an observable quantity, 
like in the case of weakly coupled Bose-Einstein condensates (30j . 

It is worth to mention that the form and corresponding frequency of a breather in an isolated chain can be obtained 
in our model only in the self-trapping breather regime, in which one can consider the limit of C — > 0. Indeed, according 
to Eq. (46), in this limit one has 8 = (z = 1) in Eqs. (19) and (20), and the breather frequency is equal to 

w(* m a», k) = u m + - 2a 2 )V 2 nax + ^ 2 max -h 2 = 2 + (|/3 - 2a)* 2 max - -k 2 . (49) 

This expression for the breather frequency is fully consistent with the known expressions for (stationary or slowly- 
moving) breather in a single a-/3-FPU chain in the small-amplitude limit, see, e.g., Refs. (3, l2l|. It is important to 
emphasize that to get this expression for breather frequency u>, one has explicitly take into account the linear increase 
in time (winding up) of the relative phase — $/2 in Eq. (19), given by Eq. (46) in the self-trapping regime. Similar 
winding up of the relative phase of two coupled macroscopic Bose-Einstein condensates in the nonlinear self-trapping 
regime has been recently measured directly in a single bosonic Josephson junction by means of interference [301 ] . This 
finding gives us an additional argument in favor of the similarity between macroscopic tunneling quantum dynamics 
and phase-coherent dynamics of weakly coupled breathers. 



III. NUMERICAL SIMULATIONS AND COMPARISON WITH ANALYTICAL PREDICTIONS 

To simulate the dynamics of coupled nonlinear lattice excitations, we numerically integrate equations (2) for two 
identical chains, ?W = ;( 2 ) = i ; qA 1 ) = qX 2 ) = a and (3^ — (3^ = (3, with absorbing edge conditions. We use the 
latter conditions in order to get rid of weak radiation, caused by the wandering breather (because such radiation 
will stay in the system forever in the case of periodic boundary conditions). For the convenience of the numerical 
simulations, the Fermi-Pasta-Ulam Hamiltonian (1) for two coupled chains with N particles can be written in the 
following equivalent dimensionless form: 



N N-l N 



,(3" 



n=l 



(50) 



where the potentials V(x) and U(x) describe, respectively, the intra- and inter-chain interactions. We normalize the 
dimensionless potential V(x) with the conditions V^O) = V'(0) = 0, V"(0) — 1, while the potential U(x) we normalize 
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with the conditions U(0) — U'(0) — 0, U"(Q) = C > 0. According to Eq. (1), we take the following form for the 
intra- and inter-chain potentials: 

V(x) = l -x 2 + l -ax 3 + X -{3x\ (51) 

U(x) = ±Cx 2 , (52) 

with C — 0.1 describing the strength of the weak inter-chain coupling. 
Then the equations of motion (2) will take the form as 

^=F{r^)-F{r ( t 1 ) + G{5^), (53) 

where r£> = - , 6$ = - , i = 1, 2, F(r) = V'(r) = r + ar 2 + fir 3 , G(S) = U'(6) = C6, 

We search for the localized nonlinear lattice excitation (discrete breather) as a solution of a set of nonlinear 

equations ^(X) = 0, where the vector X = {un\o), iin {0)} 2 ^[ n=1 gives the initial values for Eqs. (j53"|) , while the 

vector ^-"(X) = {un\t p ) — Un\o), Un\t p ) — Un\o)}^^[ n=1 gives the change of the vector X during the one period 
of breather oscillations t p — 2ir/ui. To find the value of ^(X), one needs to integrate numerically Eqs. ([53]) during 
the time interval [0, tp]. The use of such numerical method of finding of an exact breather solution is explained in 
details in Ref. 10]. The main difficulty in this numerical method is to find an appropriate initial vector X for the 
subsequent iterative solution of nonlinear equations (|53[) . 
The energy of the breather is determined correspondingly as 

2 
i=l 




A. Simulation of dynamics of discrete breathers in two weakly coupled chains 

Equations (53) for two coupled chains can be reduced to the equations for one chain in the case of symmetric, 
u n i = u n 2 = u n , and antisymmetric, u„.i = — u„,2 = u ni motion in the chains. In the symmetric case the 
Hamiltonian (50) has a usual form for a single FPU chain: 

N 1 

H = 2^{-ul + V(u n+1 -u n )}. (55) 

n=l 

For the symmetric intra-chain potential V(r), i.e. for a = 0, equations of motion in a system of two chains can also 
be reduced to the one-chain equations in the case of antisymmetric (anti-pasc) motion in the chains. In this case the 
reduced system has a standard Hamiltonian of the FPU chain on an external harmonic substrate: 

N 1 

H = 2^{~u 2 n + V(u n+1 -u n )+2Cu 2 n }. (56) 

n=l 

Discrete breathers in such chain with C > were studied in details, see, e.g., Ref. [10]. In the system with 
Hamiltonian (|56"1) . either with C = or C > 0, there are two types of discrete breathers for oj > uj max = \/4 + 2C, 
see, e.g, Refs. [3j, U, [5j. The center of symmetry of the breather coincides with a given lattice cite in one breather 
type, and it coincides with a middle point of two neighboring cites in another type. In both cases of the symmetric 
and antisymmetric breathers in the coupled chains, the vibrational energy of the nonlinear excitation, Eq. (|54p . is 
equally distributed between the chains. Numerical analysis of Eqs. (|53|) reveals that besides these (trivial) nonlinear 
eigenmodes there is an "asymmetric" mode, in which only one chain is involved in the vibration while the vibration 
of the other chain "adiabatically" follows the vibration of the main one. The characteristic form of the displacement 
patterns in one- and two-chain antisymmetric breathers is shown in Fig. 2, while the dependence of their energies on 
breather frequency to is shown in Fig. 3. 

As we can see from the comparison of Figs. 2(c) and 2(d), the magnitudes of the static displacement kinks in the 
center-of-mass displacement differences Awq CM ^ — Wq CM ^ (+oo) — w^ M \—oo) differ almost twice for the two- and 
one-chain breathers in coupled chains with cubic intra-chain anharmonic potential, in accordance with the prediction 
of Eqs. (9) and (10) that the ratio of such magnitudes is given mainly by the ratio of the total vibration energies 
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of the corresponding breathers. As follows from Fig. 3(c), the ratio of the (total vibration) energies of the two- and 
one-chain breathers for uj = 2.3, a = — 1 and (3 — 1 is indeed close to two. 

In the absence of the cubic anharmonic intra-chain potential (a = 0, (3 = 1), the energy of the antisymmetric 
two-chain breather monotonously increases from zero, for u) = Lo max = \/A + 2C ~ 2.05, to the infinity for uj — > oo, 
see Fig. 3(a). For low frequencies ui max < uj < 2.066, there is only one type of the antisymmetric two-chain breathers, 
but at ui = 2.066 = ui ij an one-chain breather mode detaches from the two-chain mode. The one-chain mode exists 
for all frequencies ui > Lu^if- As follows from Figs. 3(a) and 3(b), there is a typical bifurcation of the nonlinear mode, 
cf., e.g., Ref. [H. 

We can get an analytical estimate for the characteristic frequency u>bi / with the use of generic evolution equations 
(27) and (28). Indeed, for positive /3, <!> = ir and ui m \ = ui m 2 one has an equation for the steady-state, with 9 = 
and $ = 0, dynamical state in the coupled chains, 

C 



/2 

- max 



which has a solution for real z only for i(3^ nax > C. Using this condition and Eq. (36) for breather frequency, in the 
case of a — and k — 0, we get the following expression for uj if at the appearance of the new solution (for C = 0.1): 



uj bif = W4+-C = V425 = 2.062, (58) 

which coincides with good accuracy with the numerically revealed value u)uf — 2.066. As follows from Eq. (57), at 
the bifurcation point 3/3^ nax = C one has z = 0. This corresponds to the equal populations of the coupled chains, 
which in turn results in equal energies of the two- and one-chain breather solutions at the bifurcation, see Fig. 3 (a). 

Two types of the coupled breathers, the one- and two-chain breathers, can exist also in the presence of cubic term 
in anharmonic intra-chain potential V(x), Eq. (51). The effect of the cubic term in anharmonic intra-chain potential 
on breather dynamics in a single FPU chain was studied in Refs. @, [H, H^, [H[ • For instance, it was shown in Ref. 
that low-amplitude breathers can exist in a single FPU chain only in the case of relatively weak cubic anharmonic 
term in the intra-chain potential V(x), namely when a < -y/3/3/4. (For the considered case of (3 = 1, this requirement 
reduces to a < \/3/4 = 0.866). In fact this requirement corresponds to the requirements uj > ui m in Eq. (49) and 
Ai > in Eq. (23) for breather frequency and inverse localization length in a single FPU chain (for = and k = 0). 

From our analysis we can conclude that bifurcation of breather modes in two weakly coupled nonlinear chains with 
cubic intra-chain anharmonic potential occurs only for a < -y/3/3/8 (or a < 0.61 in the considered case). For larger 
cubic intra-chain anharmonic potential a > \/3/3/8, there are no low-amplitude breathers in the coupled chains but 
there are finite-amplitude breathers with a threshold in energy which does not go to zero for uj — > u> max + 0, see 
Fig. 3(c) and cf. with Fig. 4 in Ref. Indeed, from Eq. (23) follows that for the wandering breather, when 

(cos 6) = (sin 2 9) = 0.5, the inverse localization lengths Ai,2 are real positive only for a < y\3/3/8, when the above 
picture of the breather periodic translation (wandering) between the coupled chains takes place. On the other hand, 
for a > -y/3/3/8 (or a > 0.61 in the considered case), one has to take into account the contribution to breather 
frequency u> and inverse square localization lengths Af 2 of the finite-amplitude breather of the terms oc ^m ax (with 
positive coefficients). These terms will result in a finite (threshold) energy of the breather with minimal possible 
frequency ui m i n > ui m . According to Eqs. (14), the existence of such breather threshold amplitude and energy in the 

coupled chains with a > ^ 3/3/8 is related with the static strain in both chains {dw^ M ^ / dx — —a(\ | 2 + | ^2 | 2 )), 
which accompanies breathers in such coupled chains and which is localized around the central region of the breather, 
see Eqs. (9) and (10) and Figs. 2(c) and 2(d). 

In the coupled chains with a > -y/3/5/8, there is no bifurcation of the two-chain breather mode: there are two 
separated breather modes, the one- and two-chains ones, for all the frequencies uj > uj min , the energy of the two-chain 
breather mode is always larger than the energy of the one-chain one, and there are thresholds in energies for both 
types of breathers, see Fig. 3c. In the following we will restrict ourselves only to the case of weakly coupled /3-FPU 
chains with zero cubic anharmonic intra-chain potential and assume everywhere below that a = and (3—1. 

For the large-amplitude (high-frequency) nonlinear oscillations in two weakly coupled /3-FPU chains, there are many 
coupled (bound) antisymmetric breather modes whose strongly localized displacement patterns are shifted one with 
respect to another along the chains, see Figs. 4(a) - 4(j). In these figures the total breather energy E, Eq. (54), 
location of the breather center in the corresponding chain, and the "binding energy", i.e. the gain in energy due 
to breather coupling, AE — E[ lch ^ + E^ 1 ^ — E, are indicated for breathers with (high) frequency uj — 5. Here 
E^ 1 ^ indicate the energy of the corresponding one-chain breather located in the corresponding chain, see Figs. 4(a) 
and 4(b) for the two main types of the one-chain breathers: with the center of symmetry placed in the middle of 
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two neighboring sites or in a particular site. As it is seen from Fig. 4, the binding energy depends on the type 
of the bound breathers and apparently decreases with the increase of the inter-breather distance (along the chains). 
Besides such one- frequency coupled breather modes, there are also two- frequency nonlinear modes which bind together 
large-amplitude breather modes in two chains with different but commensurate frequencies u>\ and u>2 > tOi- The 
dependence of the energy of several types of such commensurate-frequency bound breather modes on their maximal 
frequency ui2 is shown in Fig. 5. 

We also study the interaction of breathers in two weakly coupled /3-FPU chains. For this purpose, we numerically 
integrate equations of motion (53) with the initial conditions which describe a breather with frequency u>\ on the first 
chain and a breather with frequency u>2 on the second chain. In Fig. 6 we show the change in time of the spectral 
density p(oS) of the emerged localized oscillation mode in the coupled neighboring chains with two breathers with 
initial frequencies u>\ = 2.6 and lu 2 = 3. As one can see in this figure, a two-frequency breather is generated in the 
system with the ratio of frequencies ui\ : u)2 = 6 : 7, which subsisted for 1200 units of dimensionless time (when time 
is measured in units in which ui m = V4 + C ~ 2). Later the breather with higher frequency (and energy) starts to 
adsorb the energy of the breather with smaller frequency (and energy), which decreases the frequency of the latter. In 
result, the one-chain breather with frequency to — 3.160 is formed. Hence the breather with higher energy completely 
adsorbs the breather with smaller energy. This reminds the merging of low-amplitude (small-energy) breathers in 
a single ID /3-FPU chain, which results in a single localized object (discrete breather) containing almost all of the 
initial total vibration energy of the lattice [2ll ] . If we initially prepare two breathers with more close initial frequencies 
lj\ = 2.7 and UI2 — 3, a two-chain one-frequency breather with antisymmetric displacement pattern is formed (not 
shown). In this case the interaction and energy exchange of the two breathers with initially different frequencies result 
in frequency equalization. Therefore our simulations show that despite the existence of the coupled two-frequency 
breathers, the one-frequency coupled breather modes are in general more stable. 

B. Simulation of wandering breather in two weakly coupled chains 

Here we will study the dynamics of two weakly coupled chains under excitation of a breather only in one chain 
while keeping the other chain initially at rest. This amounts in numerical integration of Eqs. (53) for chain 1 under 
(2) . (2) 

the initial conditions u n — and u n — 0. In result of inter-chain interaction, the initially located in one chain 
excitation can start to periodically translate (wander) between the coupled chains. Such inter-chain energy exchange 
can be studied quantitatively by measuring the time dependence of the energy Ei in each chain, see Eq. (54). 

In Fig. 7 we plot the energy of each chain Ei as a function of time when the immovable breather (with V = 0) is 
excited in chain 1 with a given frequency uj. For relatively small breather amplitude and corresponding dimensionless 
breather frequency, u> = 2.030 > ui m = vi + C = 2.025, an almost complete inter-chain energy exchange occurs when 
the total breather energy E = Ei + E 2 , Eq. (54), is periodically located in one of the coupled chains, see Fig. 7(a). 
From this figure we can measure the period of wandering T Wl which is close to 126. Taking into account that the 
frequency of the complete inter-chain energy exchange for small-amplitude breather is equal to C/2 = 0.05, see Eqs. 
(21) and (45), the wandering period T w is equal to Att/C = 125.66, which coincides with very high accuracy with 
the period, which can be found from Fig. 7(a). With the increase of breather amplitude and frequency, the period 
of the inter-chain energy exchange increases, in accordance with Eq. (45), but the exchange remains complete up to 
the critical breather frequency (and corresponding amplitude), which according to Figs. 7(b) and 7(c) is very close 
to lu — 2.0980. This critical breather frequency should be compared with the value of frequency of the immovable 
breather at the separatrix, u sep — \/4 + 4C = 2.0976 for C = 0.1, a = and k = 0, see Eq. (48), which also shows 
very good agreement with the prediction of the approximate physical-pendulum equation (39) (or (40)). For the 
higher breather frequencies u> > 2.098 (and corresponding breather amplitudes), the inter-chain energy exchange is no 
more complete and corresponds to the exchange of relatively small fraction of the total breather energy E, see Figs. 
7(c) and 7(d). For uj = 2.1106, the main part of the breather energy remains all the time in chain 1 were the breather 
was initially excited, which corresponds to nonlinear self-trapping of the breather, see Fig. 7(d). The transition from 
the oscillatory energy exchange, Fig. 7(b), to the self-trapped mode, Fig. 7(c), is very sharp: the change of the 
breather frequency by the crossing of the separatrix occurs in the fifth digit only. One can also see a big difference 
in energy exchange rate in the wandering breather mode, Fig. 7(a), and in the separatrix mode, Figs. 7(b) and 7(c): 
with the change of the breather frequency from u> = 2.0300 to u> = 2.098, the rate of energy exchange decreases more 
than twenty times! This observed feature of the coupled nonlinear systems is also in accordance with the prediction 
of the physical-pendulum equation (39) (or (40)) for the separatrix mode with the infinite period of the inter-chain 
energy exchange, Eq. (47). 

In Fig. 8 we show the Fourier power spectrum of the wandering breather. As it was explained above in connection 
with complex wave field of the wandering breather, Eqs. (19) and (20), the temporal Fourier spectrum of wandering 
breather consists of the main peak at breather fundamental frequency uj and its (upper and lower) frequency satellites 
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at uj ± nuJbeat, n = 1,2, where Wb ea i ~ (C — 9/3 2 \E'^ actsc / (16C))/2 is the rate of complete inter-chain energy exchange 
in the wandering breather mode, see Eq. (45). The value of LObeat ~ 0.05 for low-amplitude wandering breather in 
the chains with the coupling constant C = 0.1 is in a good agreement with the inset in Fig. 8. Therefore we can 
quantitatively reproduce both the numerically observed period, Fig. 7(a), and wandering (beating) frequency, Fig. 8, 
of the low-amplitude and macroscopically wide breather. 

The nonlinearity of the physical-pendulum Eq. (39) (or (40)) results in all higher harmonics (second, third etc.) 
of the beating frequency LObeat, and respectively in the appearance of corresponding doublets of breather frequency 
satellites: the satellites of the first and second harmonics of the beating frequency are clearly seen in the inset in Fig. 
8. As we can also see in the inset in Fig. 8, the lower-side satellites enter the continuum of low-amplitude phonons, 
u> — nuJbeat < k 'm = V 4 + C = 2.025, n — 1,2,.... The latter means that the wandering breather (weakly) emits 
low-amplitude lattice phonons, which results in (adiabatically slow) lowering of its energy. This in turn means that 
the wandering breather is not an exact nonlinear eigenmode of two weakly coupled anharmonic chains, in contrast to 
the two-chain symmetric and antisymmetric breathers. 

In Fig. 8 one can also see the third harmonic of the main breather frequency ui, which can be modelled with the use 
of the full-harmonics Eqs. (12) or (53), instead of Eqs. (14) written in the resonant approximation. For the /3-FPU 
chain (with zero cubic anharmonic intra-chain potential), there are only odd harmonics of the fundamental breather 
frequency. In a single chain, the higher harmonics of the fundamental breather frequency were first analytically 
predicted and described in Ref. 0], and were numerically observed, up to the ninth harmonic, in Ref. [22| . In 
Fig. 9 we show the power spectrum of the one-chain breather with frequency uj — 4 in the coupled chains with 
a = 0,/3 = 1, C = 0.1. This spectrum clearly shows the existence of higher odd harmonics, up to the seventh one, 
of the fundamental frequency. Since the breather with such frequency is in the self-trapped mode in the considered 
coupled system (uj = 4 > u> KK „ = 2.1), its dynamics and power spectrum are similar to the ones of the breather in a 
single chain, cf. Refs. [2l|, [22[. Importantly, the Fourier spectrum in Fig. 9 corresponds to almost "exact" breather 
solution found in the /3-FPU chain with zero vibrational background, while the breather, studied in Refs. [2~il. |22|. 
was self-created on a non-zero vibrational background in result of modulational instability of short-wavelength modes 
in the /3-FPU chain. 

In view of the present work, we can relate the observed in Ref. [22[ doublets of satellites of the fundamental breather 
frequency (and its odd harmonics) in a single chain with the beating between two breather states with different center 
locations (and symmetry), see Figs. 4(a) and 4(b). According to these figures, the breather with the inter-site center 
location, Fig. 4(a), has higher frequency than the breather with the same energy with the on-site center location, 
Fig. 4(b). According to generic Eq. (28), the frequency difference under the proper initial conditions can cause 
breather wandering (beating) between two nearest lattice sites. The beating in turn can induce the steady-state 
translation of the breather along the chain with (slow) velocity V <C 1 , proportional to the (small) beating frequency 
(as it was qualitatively explained in Ref. [22j). In Section III C we will discuss similar lateral translation of ID 
breather in a system of parallel weakly coupled nonlinear chains with transverse group velocity, proportional to the 
beating frequency in the case of two coupled chains. In contrast to the low-amplitude ID breather, laterally moving 
and losing its energy due to phonon emission in a system of nonlinear chains, see Figs. 14(a) and 14(c) below, the 
high-frequency slowly-moving breather in a single chain does not emit phonons in the chain because the fundamental 
frequency satellites (caused by the beating) are all out of the phonon band, see Fig. 3 in Ref. 22]. 

Similar picture of the complete inter-chain energy exchange takes place also for wandering breather, slowly moving 
along the chains. To model this effect, we numerically integrate Eqs. (53) with the initial conditions which correspond 
to the excitation of the moving in chain 1 breather solution, while keeping chain 2 at rest, with the following ansatz 
for lattice displacements: 

u«(i) = {-l) n V max cos(kn-ujt)/cosh[\(n-Vt)], u^(0)=0, u^(0)=0, (59) 

where ^ ma x, w = w(f mar , k) and V rs —k/uj <C 1 are the amplitude, frequency and slow velocity of the breather, cf. 
Eq. (19). 

In Figs. 10 and 11 we show the energy, versus time and site, of slowly-moving breather, with velocity V = 0.1, 
wandering between two weakly coupled chains. In Fig. 10 we show the energy of the wandering small-amplitude 
breather below the separatrix: with ^ max — 0.1, Figs. 10(a) and 10(b), and $ mra = 0.2, Figs. 10(c) and 10(d). Both 
figures show the complete inter-chain energy exchange. From Figs. 10(a) and 10(b) one can find that the period of 
the complete inter-chain energy exchange (wandering) T w « 128 for $ mnl = 0.1 is very close to the one T w 126 in 
Fig. 7: both figures describe almost harmonic wandering of a small-amplitude breather, with wandering frequency 
C/2 = 0.05. Below the separatrix, the wandering period T w increases with breather amplitude: one has T w w 179 
for ^ m ax = 0.2 in Figs. 10(c) and 10(d). Figure 11 shows the dynamics of slowly-moving breather initially excited 
in chain 1 with amplitude $ mra = 0.26, which is very close to the one at the separatrix: $m7j = ^2C/{3 = 0.2582. 
Separatrix-like dynamics, similar to the one shown in Fig. 7(b), is well established for the later delay time i>4000, 
when almost total energy of the moving breather periodically translates (wanders) between the chains. In Fig. 12 
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we show the dynamics of slowly- moving breather, initially excited in chain 1 with velocity V = 0.1 and amplitude 
^max = 0-3, which is beyond the separatrix. Here the inter-chain energy exchange is no more complete, like in the 
case of immovable breather beyond the separatrix, cf. Fig. 7(d), but the period of such exchange, T x ps 125, is much 
shorter than that close to the separatrix, cf. Fig. 11. 

Slowly-moving wandering Bose-Einstein condensate can be realized in two coupled ID atom waveguides (see, 
e.g., Refs. [ID, [EH), for the BEC of weakly- interacting atoms with non-negligible inter- waveguide tunneling coupling 
and with total initial population imbalance |z(0)| = l. 

Therefore the excitation of the low-amplitude (low-frequency) breather, either immovable or slowly-moving, in 
one chain always results in its periodic transverse translation (wandering) between the two weakly coupled chains. 
Hence a natural question arises what happens if one deals with a system of M > 2 parallel coupled anharmonic 
chains. Whether the wandering of the breather, initially excited at the edge (outermost) chain, across all the chains 
is possible? To study this problem, we start with a system of M > 2 parallel weakly coupled anharmonic chains. 



C. Simulation of breathers in a system of M > 2 parallel weakly coupled anharmonic chains 

Quasi- ID system of M parallel weakly coupled anharmonic chains, with nearest-neighbor intra- and inter-chain 
interactions, is described by the following Hamiltonian: 

M N M N—l M-l N 

H = E E 5"m,n + E E y ( ) + E E U(u m+ x t „ - Um >n ), (60) 

m— 1 n—l m— 1 n—l rn — 1 n—l 

where V(x) and U(x) are given by Eqs. (51) and (52), and n = 1, N and m = 1, M numerate, respectively, sites 
along the chains and the chains. 

Hamiltonian (60) generates corresponding equations of motion, 



dH 



(61) 



OUrn.n 

which in the linear approximation have the form: 

^m,n — ^m,n+l 2ti m n -)- W m , n _ 1 -f- C(u m ^i jn 2u mn -)- llm— l,n)' (^^) 

Plane linear waves (phonons) in such system, with 

Um,n — uexp[iqin + %qi,m — iut], (63) 

have the dispersion: 



^(91,92) = V 2 ! 1 - cos 9i + c ( l - cos 92)], (64) 

where both the intra-chain lattice period and inter-chain spacing are taken equal to unit, and therefore Q < (gi, q%) < tt. 
Minimal phonon frequency in this translationally-invariant system is zero, w(0, 0) = 0, while cut off phonon frequency 
is w(7T, 7r) = 2V1 + C (which is equal to u(tt, it) = 2.0976 » 2.1 for C = 0.1). 

Since we will be interested in the short-wavelength excitations, with qi w n, the corresponding phonon frequency 
for C < 1, 

u(q 1 ^7r,q 2 ) = 2 + Csm 2 (^q 2 ), (65) 
determines the phonon group velocity across the chains: 

Vj.(gi w 7r,g 2 ) = 5 = -Csm 92 • (66) 

dg 2 2 

Now we turn to the nonlinear dynamics of M weakly coupled parallel nonlinear chains with Hamiltonian (60). 
We will integrate Eqs. (61) with the initial condition which describes exact discrete breather in the m-th chain 
(1 < m < M ) under the condition of immovability of the rest of the chains, to study the time dependence (for t > 0) 
of the vibration energy in the chains: 



^m,n ~l~ Vrn.n + Vm,n—\ "b U m , n + U m — l,n) — ^ ' E n m , (67) 



71=1 
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where V m . n = V(u m . n+1 -u m . n ), U m<n = U (u ro +i, n -ii m ,„), and breather energy E = J2m=i E ™ = E n =i Em=i ^n,m- 
By measuring the time dependence of E n m , we can define the breather center location (n c ,m c ), 

M N M N 

n c = ^ ^ n Pn,m, m c = ^2 X! m Pn,m, (68) 

m—1 n—1 m — 1 n—1 

where p n ,m — E n>m /E describes the distribution of normalized breather energy in 2D lattice. With these breather 
quantities we can measure the longitudinal D x (along the chain axis) and transverse D y breather diameters: 



M N 
m—1 n—1 



1/2 



Dy 



M N 

- m c ) 2 p n ^ 

m—1 n—1 



1/2 



(69) 



In Figs. 13 (a,b,c,d) we show the time dependence of energies of the first, i = 1, and the last, i = M, coupled 
chains for M = 2, 3, 4, 5, for the time interval just after breather excitation in the first chain, left panel, and for the 
later time, right panel. In the case of M = 2, there is a periodic (harmonic) complete energy exchange between the 
first and the second chains, see Fig. 13(a). In the case of M = 3, there is periodic (non-harmonic) recurrence of the 
complete energy accumulation in the first chain, and the period of such recovery is twice larger (the recurrence rate is 
twice smaller) than that in the case of M = 2. For M = 3, the time dependence of the first chain energy recurrence 
can be roughly approximated as cos 4 (Ci/8) (instead of cos 2 = cos 2 (Ci/4) for M = 2, Eqs. (19) and (45)), similar 
to the time dependence of the population of the initially excited state in a three-state (spin 1) atomic system, see, 
e.g., Ref. [56| . Here C/4 plays role of the Rabi frequency, which is twice smaller than the rate of the complete energy 
exchange (energy recurrence) in the case of M = 2, when it is equal to (7/2, see Eq. (45). In the case of M = 4 
and M = 5, the recurrence of energy of the first and the last chains becomes quasi-periodic, see Figs. 13 (c) and 13 
(d), but still the (approximate) period Tm of such recurrence scales with the number of chains M as Tm oc (M — 1). 
For M > 6, C = 0.1, the initially localized (in chain 1) excitation spreads its energy in the whole system of weakly 
coupled chains. 

The dependence of the recurrence period Tm on M we can relate with the transverse group velocity, Eq. (66), of the 
high-frequency phonons (with q\ rs 7r) in the system of weakly coupled anharmonic chains. The transverse wavevector 
<72 in Eq. (66) is equivalent to the relative phase of the neighboring chains. In the regime of almost-harmonic energy 
transfer between the neighboring chains, the relative phase is always close to tt/2, see Fig. 1. It means that the 
transverse wavevector qi in the expression (66) for the transverse group velocity of the wandering breather should 
also be (approximately) equal to tt/2. In the case of 2D system of coupled oscillator chains, q^ — tt/2 corresponds to 
the case when, say, only the odd chains are excited at a given moment while their nearest neighbors, even chains, are 
at rest. This gives V± « |C for the transverse group velocity. Importantly, this characteristic group velocity does 
not depend on the number of the coupled chains M. Therefore we can estimate the period of the first chain energy 
recovery as Tm — 2A(M — 1)/Vj_ ~ &A(M — 1)/C with some dimensionless factor A, which is consistent with our 
numerical observation (with Af»3). 

The same transverse phonon group velocity one can estimate from Fig. 14 as the speed of breather spreading 
across the chains. Figure 14 shows the time dependence of energy distribution among the chains when the breather is 
initially excited in the edge chain, Figs. 14(a) and 14(b), or in the central chain, Figs. 14(c) and 14(d), in the system 
of M = 50 coupled chains with N = 50. As follows from Figures 14(a) and 14(c), the initial breather energy has 
spread for 20 chains for approximately 500 time units. This gives us quantitative estimate of 0.04 for the transverse 
group velocity, which is rather close to our analytical estimate with the use of Eq. (66) for q^ rs tt/2: V± ~ \C = 0.05. 
Figure 14 also shows that the appearance of localized breather in a system of coupled chains, with M ^> 2, has a 
threshold in breather frequency, tothresh — 2.15 in Fig. 14(b) and oothresh = 2.17 in Fig. 14(d), similar to the case of 
two coupled chains. The threshold breather frequency, which corresponds to the appearance of localized breather in 
a system of coupled chains, we should compare with breather frequency at the separatrix in a system of two coupled 
chains, uj sep w 2.1, see Eq. (48) for C = 0.1, a — 0, k = and Figs. 7(b), 7(c) and 11. Since all the above 
frequencies are close, we can get the estimate for the threshold breather amplitude for its localization in 2D system 
of weakly coupled chains: ^ r a e x sh - ^fC/P- Fot ^max < *mS sh > the 1D breather, which was initially excited in one 
chain, will start to translate laterally to the neighboring chains, will spread its energy among them and lose it due to 
phonon emission, and finally will decay into small-amplitude phonons due to lowering of its frequency up to the cut 
off phonon frequency io ma x — 2\J\ + C, see Figs. 14(a) and 14(c). Such evolution of low-amplitude ID breathers can 
also be related with the above-mentioned conclusion that wandering breather is not an exact solution of the nonlinear 
system even in the case of two coupled anharmonic chains. In contrast to this behavior of low-amplitude breathers, 
ID breather with the amplitude ^ ma x > ^mLx 3h i s self-trapped and remains localized mainly in the chain of its initial 
excitation, see Figs. 14(b) and 14(d). [In Fig. 14 (d) one can also see a partial (incomplete) energy exchange between 
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the central chain and its nearest neighboring chains, similar to the incomplete energy exchange in two coupled chains 
beyond the separatrix, cf. Figs. 7 (d) and 12]. This phenomenon resembles the so-called delocalizing transition 
in 2D system, when the wave field abruptly changes its character from spatially localized to the extended one, cf. 
similar delocalizing transition for polarons in 2D and 3D lattices [35[ and for Bose-Einstein condensate in 2D optical 
lattice [361 ] . In our case, the delocalizing transition occurs by the decrease of the initial breather amplitude $ mra (or 
frequency lS) from the value $> max > ^ r a % sh ~ yJC/fi (or ui > uj thresh ) to the value 9 max < ^ r a % sh (or w < u) thresh ). 
Such transition is related with finite energy threshold for the creation of solitons and breathers in 2D and 3D systems, 
see Refs. |57j,|58|, and is absent in ID (/3-FPU or discrete nonlinear Schrodinger equation [35|,l36j) systems. Indeed, the 
threshold breather amplitude \f r ^£J s ' 1 m strongly anisotropic quasi-lD system vanishes in the limit C — * of a single 
ID chain since ^max Sh ~ V & / (3- Similar threshold breather amplitude ^^ms ~ \/C//3 for breather derealization 
should also appear in 3D array of parallel weakly coupled nonlinear chains, with coupling constant C -C 1, which are 
described by the FPU Hamiltonian, similar to the one given by Eq. (1). 

To get more deep insight into the delocalizing transition of discrete breather in 2D system of weakly coupled chains, 
we confirm numerically that discrete breather in 2D system of weakly coupled chains is stable only in the frequency 
range [ub, oo), where u>b = 2.12 for C = 0.1 is higher than the cut off phonon frequency lj(tt, it) = 2yl + C « 2.1. This 
means that there is frequency gap [u)(tt, 7t), cjg] above the cut off phonon frequency in which there are no solutions for 
localized breathers in 2D system: breathers with frequencies in the gap spread across all the 2D lattice. The frequency 
u>b is close to the aforementioned and indicated in Fig. 14 threshold frequencies uJthresh-, but the frequency ujf, does 
not depend on the way of the breather excitation in 2D lattice, in contrast to the threshold frequencies. In Fig. 15 we 
plot the breather energy distributions along and across the chains for the one-chain, Fig. 15(a), and antisymmetric 
(anti-phase) two-chain, Fig. 15(b), breathers in 2D lattice with to = ujb = 2.12, which show that transverse localization 
length of the discrete breather has the order of the inter-chain spacing for u> = u>b (and further decreases for u> > u>b, 
see Fig. 16 (b) below). 

In addition to the one- and two-chain breathers, there are also three-chain, four-chain etc. in- and anti-phase 
breathers in a multi-chain system. In Fig. 16 we show the vibration breather energy and the measured (with the use 
of Eqs. (67)-(69)) longitudinal and transverse breather diameters for the one-, two-, three- and four-chain anti-phase 
breathers. As one can see from this figure, breather energy grows as uo A in the high-energy limit. This dependence 
can easily be understood from Eqs. (16) and (22) for breather energy Eb and frequency u: one has Eb oc ^m ax and 
u) oc >t m(a , in the large-amplitude limit, see also Ref. 22]. Breather with threshold frequency ui = uib has minimal 
energy: Eb = 0.517 and Eb = 0.740 for the one- and anti-phase two-chain breathers, respectively. From Fig. 16 (b) 
we can also see that, in accordance with Fig. 15 (a), the transverse localization length of the one-chain breather at 
the threshold frequency is very close to the inter-chain spacing and further decreases with the increase of breather 
energy. [Both longitudinal and transverse breather diameters saturate at the amplitude-independent values in the 
high-energy limit, see, e.g., Q, when the transverse diameter of the m-chain breather simply is D y « B(m + 1) 
with B = 1, see Fig. 16(b).] This means that there are no laterally-localized envelope-soliton vibrational excitations 
in our 2D system of parallel coupled nonlinear chains, in contrast to the laterally-localized envelope-soliton optical 
excitations in 2D arrays of parallel coupled nonlinear optical waveguides [TJ]. Therefore the inter-chain breather 
translation and wandering are possible only in a system of small number of coupled chains (M < 5 for C — 0.1), while 
in a multi-chain 2D system (M > 6 for C = 0.1) the vibrational breathers can move only along the chains. Such ID 
motion of ID breathers in 2D system of multiple parallel weakly coupled oscillator chains resembles the ID motion 
of ID Bo se g ases in 2D array of multiple parallel weakly coupled atom waveguides (a quantum Newton's cradle) [28[ 
(see also |59j). 

IV. SUMMARY 

In summary, we have found, both analytically and numerically, two qualitatively different regimes of energy ex- 
change between phase-coherent breathers (intrinsically localized short-wave nonlinear excitations) in two weakly linked 
nonlinear chains. In the low-amplitude mode, the breather performs periodic transverse translations (wandering) be- 
tween the coupled chains. In the large- amplitude mode, the breather is self-trapped in one chain. These two breather 
modes are detached by a separatrix, at which the rate of the inter-chain energy exchange vanishes. These regimes 
have a profound analogy, and are described by a similar pair of equations, to the Josephson plasma oscillations and 
nonlinear self-trapping, recently observed in a single bosonic Josephson junction [3fj| . The predicted evolution of the 
relative phase of two weakly coupled coherent breathers in both regimes is also analogous to the evolution of relative 
quantum mechanical phase between two weakly coupled macroscopic condensates, which was directly measured in a 
single bosonic Josephson junction by means of interference [3(j. On the basis of this profound analogy, we predict a 
new tunneling regime of two weakly linked Bose-Einstein condensates in which their relative phase oscillates around 
^ modulo 7r, which can be observed by means of interference. The similarity between the classical phase-coherent 
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excitation exchange and macroscopic tunneling quantum dynamics found here can encourage new experiments in both 
fields. 

We also show that static displacements of the coupled chains with nonlinear localized excitation, induced by cubic 
term in the intra-chain anharmonic potential, scales approximately as the total vibrational energy of the excitation, 
either one- or two-chain one, and do not depend on the inter-chain coupling. This observation is also valid for a narrow 
stripe of several parallel coupled nonlinear chains. We also study two-chain breathers, which can be considered as 
bound states of discrete breathers with different symmetry and center locations in the coupled chains, and bifurcation 
of the anti-phase two-chain breather into the one-chain one. Bound states of two breathers with different commensurate 
frequencies are found in the two-chain system. Merging of two breathers with different frequencies in one breather 
in two coupled chains is observed. Wandering of low-amplitude breather in a system of several, up to five, coupled 
nonlinear chains is studied, and the dependence of the wandering period on the number of chains is analytically 
estimated and compared with numerical results. 

Dclocalizing transition of ID breather in 2D system of a large number of parallel coupled nonlinear chains is 
described, in which the breather, initially excited in a given chain, abruptly spreads its vibration energy in the 
whole 2D system upon decreasing breather frequency or amplitude below the threshold one. The threshold breather 
frequency is above the cut off phonon frequency in 2D system, the threshold breather amplitude scales as square root of 
the inter-chain coupling constant, and breather vibration energy is localized mainly in one chain at the derealization 
threshold. Similar delocalizing transition for ID breather should also occur in 3D array of parallel coupled nonlinear 
chains. Dclocalizing transition of discrete vibrational breather in 2D and 3D systems of coupled nonlinear chains has 
an analogy with delocalizing transition for Bose-Einstcin condensates in 2D and 3D optical lattices. 
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FIG. 1: Phase portrait of wandering breather in two weakly coupled nonlinear chains or two weakly linked BECs in a symmetric 
double-well potential with the initial conditions |z(0)|=l and $(0)=7r/2 modulo tv, which is given by (kcgs"! 1 ) 2 + z 2 — 1, 
k = 2C/3/3't> 2 nax , Eq. (35). Lines 1 - 8 correspond, respectively, to k = 0.5, 0.8, 0.95, 1, 1.05, 1.25, 2, and 10. Lines 1 -3 
describe the self-trapped mode, line 4 describes the separatrix, lines 5-8 describe the wandering breather or new tunneling 
mode of two weakly linked BECs. 




FIG. 2: (color online) Displacement patterns for (a) one-chain breather and (b) two-chain antisymmetric breather with frequency 
lo — 2.3 for the quartic intra-chain anharmonic potential (a = 0, (3 = 1); (c) and (d) show displacement patterns for one-chain 
and two-chain breathers, respectively, in the presence of cubic and quartic intra-chain anharmonic potentials (a = — 1, j3 = 1) 
for the same breather frequency. The patterns are shown at the moments when all particle velocities are zero. 
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FIG. 3: (color online) Dependence of breather energy E in a system of two coupled chains on breather frequency u for different 
strengths of cubic intra-chain anharmonic potential: a — (a), a — —0.5 (b) and a = — 1 (c), with the same strength of quartic 
anharmonic potential /3 = 1. Lines 1, 3, 5 correspond to antisymmetric two-chain breather, lines 2,4, 6 correspond to one-chain 
breather. All the lines correspond to the breathers with centers of symmetry located in the middle of two nearest particles, 
cf. Figs. 2(a) - 2(d). The dotted vertical line in all the plots shows the upper boundary of phonon spectrum in the two-chain 
system 74 + 2C » 2.05. 
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FIG. 4: (color online) Displacement patterns for two-chain breathers with u — 5 in a system of two chains with pure quartic 
intra-chain anharmonic potetial, with a = 0, /3 = 1. The patterns are shown at the moments when all particle velocities are 
zero. Here E is breather energy, n\ and n\ are the locations of the centre of symmetry of the breather in the i-th chian, 

AE = E (lch) + E (lch) _ Eis 

gain in energy due to breather coupling (breather "binding energy"). 
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FIG. 5: (color online) Breather energy E in a system of two chains versus frequency u>. Dashed line 1 in the inset corresponds 
to symmetric two-chain breather, line 2 describes two-chain antisymmetric breather, line 3 describes one-chain breather. Lines 
4 - 7 in the main plot correspond to two-chain excitations which couple breathers with two different commensurate frequencies 
with ratios u)\ : 102 — 1 : 2, Wj : L02 = 2 : 3, u)\ : L02 = 3 : 4 wi : L02 — 4 : 5, respectively. Frequency to in the plot always denotes 

UJ2- 
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FIG. 6: Time evolution of spectral density p(w) of nonlinear localized excitation of two weakly coupled chains under the initial 
excitation of a breather with frequency oj\ — 2.6 in one chain and breather with frequency L02 = 3 in another chain. 
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FIG. 7: (color online) Time evolution of chain energies E\ (red) and E2 (blue) of immovable breather in chains 1 and 2 versus 
time t, obtained from numerical solution of Eqs. (53) with the initial breather excitation in chain 1 (for immovable chain 2) 
with frequency (a) w = 2.0300, (b) lo = 2.0980, (c) w = 2.0981, and (d) lo = 2.1106. Two identical chains with a = 0, /3 = 1, 
C = 0.1 and absorbing edges were used in simulations. Separatrix solution of Eq. (39) corresponds to u = 2.0976. Time 
evolution is shown after t = 10 5 from the excitation instant when all transients died out. 
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FIG. 8: Spectral density p(u>) of wandering breather in a system of two weakly coupled chains under the initial excitation of a 
breather with frequency cu = 2.05. Dotted vertical line indicates the upper bound of phonon frequency cu m = V^TC ^ 2.025 
in Eqs. (14). Inset shows fine details of the spectral density in close vicinity of the fundamental breather frequency. 
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FIG. 9: Spectral density p(ui) of one-chain large-amplitude breather with frequency u> — 4 (in self-trapped regime), in two 
coupled chains with C = 0.1, a = 0, P = 1. 
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FIG. 10: Periodic translation ol slowly-moving breather between two coupled chains with C = 0.1,a = 0,/3=1. Time 
dependence of breather energy distribution in chain 1, (a) and (c), and chain 2, (b) and (d), is shown. The breather was 
initially excited in chain 1 with immovable chain 2, with velocity V = 0.1 and amplitude ^ ma x = 0.1, (a) and (b), and with 
velocity V — 0.1 and amplitude ^ ma x = 0.2, (c) and (d). 
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FIG. 11: Energy of slowly-moving wandering breather close to the separatrix in chain 1 (a) and chain 2 (b) versus time and 
site. The breather is initially excited in chain 1 with immovable chain 2, with velocity 1^=0.1 and amplitude ty ma x=0.26. The 
separatrix solution of Eq. (|40|l corresponds to ^' mal =0.2582. Separatrix-like dynamics, similar to the one shown in Fig. 7(b), 
is well established for t>4000. 
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FIG. 12: Distribution of breather energy in chain 1 (a) and chain 2 (b) versus time in the self-trapped regime of slowly-moving 
breather in two coupled chains with C = 0.1, a = 0, (3 = 1. The breather was initially excited in chain 1 with immovable chain 
2, with velocity V = 0.1 and amplitude ^ ma x = 0.3. 
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FIG. 13: (color online) Time dependence of breather energy in the first (red) and last (blue) chains in a system of 2 (a), 3 (b), 
4 (c), and 5 (d) chains with C = 0.1, a = 0, /3 = 1. The breather is always initially excited with frequency to — 2.05 in chain 
1, with immovable rest of the chains. 




FIG. 14: Time dependence of breather energy distribution between chains in a system of 50 coupled chains with C = 0.1, 
a = 0, /? = 1. Here m is chain number and E m is energy of the m-th chain. The breather was initially excited with frequency 
u) = 2.14 in chain 1 (a), with frequency u> — 2.15 in chain 1 (b), with frequency u> — 2.16 in chain 25 (c), with frequency 
u) — 2.17 in chain 25 (d), with immovable rest of chains. 
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FIG. 15: Distribution of vibration energy in 2D lattice of weakly coupled nonlinear chains, with a = 0, (3 = 1 and C = 
along (n) and across (m) the chains for one-chain (a) and anti-phase two-chain (b) breather with frequency lo = 2.12. 
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FIG. 16: Dependence of breather vibrational energy (a), and longitudinal D x and transverse D y breather diameters (b) versus 
breather frequency for one-chain (curves 1,5,6), anti-phase two-chain (curves 2,5,7), anti-phase three-chain (curves 3,5,8), and 
anti-phase four-chain (curves 4,5,9) breathers in a multi-chain system (M = 20, -/V = 50). Dashed line corresponds to the 
longitudinal breather diameter. Vertical dotted line shows the cut off phonon frequency. 



